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INTRODUCTION 


To  increase  the  maximum  elastic  carrying  capacity  and  to  enhance  the 
fatigue  life,  residual  stresses  are  often  produced  in  tubes  through  autofrettage 
(ref  1).  Many  solutions  are  reported  for  the  hydraulic  autofrettage  process. 

The  thick-walled  cylinders  are  subjected  to  uniform  internal  pressure  of  suf¬ 
ficient  magnitude  to  cause  plastic  deformation  and  then  the  pressure  is  removed. 

A  more  economical  way  of  producing  residual  stresses  in  thick-walled  cylin¬ 
ders  is  the  swage  autofrettage  process.  This  process  is  carried  out  by  a 
mandrel,  the  diameter  of  which  is  greater  than  the  inner  diameter  of  the  tube. 
The  mandrel  is  driven  through  the  tube  from  one  end  to  the  other.  A  rigorous 
analysis  of  this  process  is  difficult.  However,  recently  a  simple  analysis  of 
the  swage  autofrettage  process  was  reported  (ref  2).  The  model  used  was  a 
simplified  plane-strain  problem  of  the  mandrel-tube  assembly.  The  steel  tube 
was  assumed  to  be  elastically-ideally  plastic,  obeying  Tresca's  yield  criterion 
and  the  associated  flow  theory,  but  the  tungsten  carbide  mandrel  was  elastic. 

The  deformation  and  stress  distribution  during  swaging  were  obtained  by  solving 
the  shrink-fit  problem  beyond  the  elastic  limit.  After  swaging,  the  permanent 
bore  enlargement  and  residual  stresses  were  calculated  by  an  unloading  analysis 
(ref  2),  taking  into  account  the  8auschinger  effect  and  the  strain-hardening 
during  unloading  (ref  3). 

The  solution  reported  in  Reference  2  is  in  closed-form  and  the  numerical 
results  indicate  that  the  agreement  between  the  calculated  and  experimental  data 
is  excellent  in  zones  with  larger  wall  ratios  but  not  so  good  in  zones  with  wall 
ratios  less  than  two.  The  differences  in  thinner  sections  may  be  due  to  the 

References  are  listed  at  the  end  of  this  report. 
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longitudinal  bending  effect  since  the  simplified  analytical  analysis  is  one¬ 
dimensional  and  bending  is  neglected.  In  order  to  determine  the  longitudinal 
effect,  a  two-dimensional  analysis  based  on  the  finite  element  method  has  been 
conducted.  In  this  report,  the  finite  element  solutions  are  presented  for  both 
one-  and  two-dimensional  models  and  a  comparison  of  the  results  is  given. 

METHOD  OF  ANALYSIS 

Since  the  total  length  of  che  tube  is  about  sixty  times  the  diameter  of  the 
mandrel,  a  complete  finite  element  analysis  of  the  swage  autofrettage  process  is 
very  difficult.  As  the  mandrel  is  driven  through  the  tube  from  one  end  to  the 
other,  the  simulation  requires  the  study  of  elastic-plastic  moving  contact  and 
separation  history  between  two  deforming  bodies.  In  addition,  considerable 
amounts  of  computer  storage  and  run  time  are  required.  In  the  present  study, 
however,  approximate  finite  element  models  were  chosen  to  represent  swaging  in 
only  a  part  of  the  tube  (zone  3).  We  consider  the  process  as  quasi-static  and 
neglect  the  effect  of  sliding  and  friction  between  the  mandrel  and  tube.  We 
want  to  obtain  the  information  about  the  deformations  and  stresses  for  a  section 
at  only  two  particular  stages,  i.e.,  when  the  mandrel  is  at  or  far  away  from  the 
position  of  interest.  To  achieve  this  purpose,  we  can  simplify  the  simulation 
by  studying  two  related  problems,  i.e.,  shrink-fit  and  complete  unloading. 

When  the  mandrel  is  at  the  position  of  interest,  we  consider  a  shrink-fit 
problem  of  the  mandrel-tube  assembly  to  obtain  the  maximum  deformation  and 
stresses  during  swaging.  When  the  mandrel  is  driven  far  away  from  the  section, 
we  study  it  as  a  complete  unloading  problem  of  the  mandrel-tube  assembly  to 
obtain  the  information  about  the  permanent  bore  enlargement  and  residual 
stresses  after  swaging.  Figure  1  shows  a  one-dimensional  interference-fit 
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problem  of  the  mandrel-tube  assembly.  Initially,  the  inner  and  outer  radii  of 
the  tube  are  a  and  b  and  the  radius  of  the  mandrel  is  c.  Given  the  interference 
I  »  c-a,  we  can  determine  the  interference  pressure  p  and  the  deformation  and 
stresses  in  the  mandrel  and  tube.  In  general,  this  problem  can  be  solved  only 
by  an  iterative  approach.  If  the  mandrel  were  rigid,  then  the  direct  approach 
using  displacement  constraints  could  be  applied.  The  results  based  on  this 
approach  were  obtained  so  we  could  discuss  the  effect  of  elasticity  in  the 
mandrel.  The  actual  strength  ratio  of  tungsten  carbide  to  steel  is  about  three. 
For  the  problem  considered  here,  it  is  reasonable  to  assume  that  the  steel  tube 
is  elastic-plastic,  obeying  Mises'  yield  criterion  and  the  associated  flow 
theory,  but  the  tungsten  carbide  mandrel  remains  elastic.  The  finite  element 
analysis  is  carried  out  by  using  the  nonlinear  program,  ABAQUS  (ref  4).  Two 
types  of  elements  used  are  shown  in  Figure  2.  The  axi symmetric  solid  elements 
( CAX4 )  are  used  to  model  the  tube  and  mandrel.  The  interface  elements  ( INTER2A) 
are  used  to  model  the  separation  or  interference  fit  between  the  mandrel  and 
tube.  Truss  elements  (CID2)  can  also  be  used  to  model  the  mandrel  because  the 
displacement  Uj  is  directly  related  to  the  external  pressure  p  by 

Ui/c  =  -(l-v1-2v1*)p/E 

where  Ej,  vj  are  elastic  constants  of  the  mandrel. 

FINITE  ELEMENT  MODEL 

Figure  3  shows  a  two-dimensional  finite  element  model  (E3)  chosen  to  repre¬ 
sent  tne  swaging  process  in  zone  3.  The  model  is  considered  symmetric  with 
respect  to  z  =  0  so  that  only  half  of  the  model  is  shown.  We  have  used  133  and 
21  elements  of  type  CAX4  to  represent  the  tube  and  mandrel,  respectively,  with 
a  =  1,  b  =  1.431,  c  =  1.007415.  There  are  eight  interface  elements  of  type 
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INTER2A  to  represent  the  interaction  between  the  tube  and  mandrel.  Figure  3a 
shows  an  interference-f it  problem  of  the  mandrel-tube  assembly.  This  model  is 
used  to  determine  the  maximum  deformation  and  stresses  during  swaging.  Figure 
3b  shows  a  complete  unloading  problem  when  the  two  parts  are  separated.  This 
problem  is  used  to  determine  the  permanent  deformation  and  residual  stresses 
after  swaging.  In  order  to  determine  the  longitudinal  bending  effect,  we  would 
like  to  compare  the  two-dimensional  analysis  with  the  one-dimensional  analysis. 
The  one-dimensional  model  (El)  consists  of  ten  elements  (of  type  CAX4)  each  for 
the  tube  and  mandrel  with  one  interface  element  (of  type  INTER2A).  Another  one¬ 
dimensional  model  (E2)  has  the  mandrel  represented  by  one  or  two  truss  elements 
of  type  CID2.  The  material  constants  used  are  E  =  200,  a0  =  1,  v  =  0.3  for  the 
high  strength  steel  and  Ej  =  590,  =  3.33,  =  0.258  for  the  tungsten  car¬ 

bide.  The  materials  exhibit  no  strain-hardening.  In  the  modeling  and  com¬ 
putation  we  have  used  the  dimensionless  quantities  with  the  inner  radius  (2.283 
inches)  as  the  unit  length  and  the  initial  yield  stress  (150  Ksi)  as  the  unit 
stress.  The  actual  quantities  can  be  obtained  easily  if  needed. 

In  the  above  three  models  (El,  E2,  E3),  the  tube  is  elastic-plastic,  but 
the  mandrel  remains  elastic.  If  the  strength  ratio  of  the  mandrel  material  to 
tube  material  is  very  large,  then  the  mandrel  can  be  regarded  as  rigid.  In 
order  to  determine  the  effect  of  elasticity  in  the  mandrel,  we  have  chosen  three 
finite  element  models  (Rl,  R2,  R3).  Models  R1  and  R2  represent  one-dimensional 
plane-strain  and  plane-stress  cases,  respectively.  We  have  used  ten  elements  of 
type  CAX4  to  represent  the  tube.  The  model  R3  is  the  same  as  the  model  E3  shown 
in  Figure  3  except  that  the  mandrel  is  replaced  by  a  rigid  block. 

Following  the  instructions  given  in  Reference  4,  we  have  prepared  the  input 
data  for  each  of  the  six  finite  element  models.  For  each  model  we  ran  the 


problem  in  two  steps,  i.e.,  loading  and  unloading.  The  input  deck  for  the 
finite  element  analysis  of  model  El  is  shown  in  Table  I. 

RESULTS  AND  DISCUSSIONS 

For  each  of  the  six  models  (Rl,  R2,  R3,  El,  E2,  E3)  discussed  in  the  pre¬ 
ceding  section,  we  have  run  the  finite  element  program  successfully.  The 
numerical  results  for  the  displacements,  strains,  and  stresses  in  the  tube 
during  and  after  swaging  have  been  obtained.  Only  the  results  for  the  stresses 
along  the  radial  direction  near  z  =  0  and  the  displacements  along  the  bore  are 
presented  graphically. 

When  the  mandrel  is  assumed  to  be  rigid,  the  displacement  at  the  bore  is 
equal  to  the  given  interference.  The  results  for  the  stresses  based  on  models 
(Rl,  R2,  R3)  are  presented  in  Figures  4  through  6.  When  the  interference  is 
only  half  of  the  maximum,  the  state  of  stresses  remains  elastic  as  shown  in 
Figures  4  and  5.  When  the  maximum  interference  (I  =  0.007415}  is  reached,  the 
state  of  stresses  is  elastic-plastic.  The  effect  of  interference  on  the  distri¬ 
butions  of  hoop  and  axial  stresses  can  be  seen  in  Figures  4  and  5,  respectively. 
By  comparing  the  results  for  model  Rl  (one-dimensional,  plane-strain  case)  and 
model  R3  (two-dimensional  case),  we  can  also  see  the  influence  of  the  longitudi¬ 
nal  effect  on  the  hoop  and  axial  stresses.  The  influence  on  the  maximum  axial 
stresses  is  very  significant  as  shown  in  Figure  5.  Unloading  after  the  maximum 
interference  is  reached,  we  have  obtained  the  residual  stresses  as  shown  in 
Figures  5  and  6  for  the  axial  and  hoop  stresses.  A  comparison  of  these  residual 
stresses  indicates  that  the  differences  between  one-  and  two-dimensional  models 
(Rl  and  R3)  are  very  minor.  Models  Rl  and  R2  represent  plane-strain  and  plane- 
stress  cases,  respectively,  and  both  models  are  one-dimensional. 


When  the  mandrel  is  considered  as  elastic,  the  interference-fit  assembly  is 
solved  iteratively.  The  same  results  for  the  one-dimensional  models  (El  and  E2) 
have  been  obtained.  A  comparison  of  two  models  (El  and  Rl)  for  the  hoop 
stresses  during  and  after  swaging  is  shown  in  Figure  7.  The  elasticity  in  the 
mandrel  reduces  the  amount  of  overstrain  from  70  to  60  percent.  The  numerical 
results  for  the  two-dimensional  model  (E3)  are  presented  in  Figures  8  through 
11.  Figure  8  shows  the  distributions  of  hoop  stresses  during  and  after  swaging. 
Figure  9  shows  the  corresponding  distributions  of  maximum  and  residual  axial 
stresses.  Also  shown  in  Figures  8  and  9  are  the  one-dimensional  results  based 
on  model  El.  A  comparison  of  the  results  based  on  models  El  and  E3  can  deter¬ 
mine  the  two-dimensional  effect  on  these  stresses.  In  Figure  10  we  show  the 
results  for  the  radial  stresses  based  on  four  models  (El,  E3,  Rl,  R3).  Finally, 
the  results  based  on  several  models  for  the  radial  displacement  along  the  bore 
are  presented  in  Figure  11.  The  displacements  during  and  after  swaging  are 
represented  by  U  and  U" ,  respectively.  Also  shown  in  the  figure  is  the  measured 
permanent  bore  enlargement.  By  comparing  the  results  based  on  models  El  and  Rl, 
the  elasticity  effect  gives  a  smaller  value  for  U".  If  we  include  the  two- 
dimensional  effect  with  model  E3,  we  get  a  value  for  U"  even  smaller  than  that 
based  on  the  one-dimensional  model. 
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TABLE  I.  THE  FINITE  ELEMENT  INPUT  DECK  FOR  MODEL  El 


★HEADING 

TUBE-MANDREL  ASSEMBLY  AND  SEPARATION 
★NODE 

1.. 

II, 1.007415 

21,1.0 

31.1.431 

101..  ,  0.05 

III, 1.007415,  0.05 

121,1.0  ,  0.05 

131.1.431  ,  0.05 

★NGEN , NSET =SIDE1 
1,11 

101,111 

*NGEN,NSET=SIDE2 

21,31 

121,131 

★NSET , NSET =BORE 
1,101 

★ELEMENT, TYPE=CAX4 
1,1,2,102,101 
11,21,22,122,121 
★ELGEN , ELSET=MANDREL 
11,10 

★ELGEN, ELSET=TUBE 

11,10 

★SOLID  SECTION , ELSET=MANDREL , MATERIAL=CARBIDE 

★MATERIAL , NAME=CARBIDE 

★ELASTIC 

590.,  .258 
★PLASTIC 
3.33 

★SOLID  SECTION, ELSET=TUBE ,MATER I AL=STEEL 
★MATERIAL , NAME=STEEL 
★ELASTIC 
2.E2,  .3 
★PLASTIC 
1. 

★ELEMENT, TYPE=INTER2A,ELSET=SFIT 
101,111,11,121,21 
★INTERFACE , ELSET=SFIT 
★FRICTION 
.0 

★BOUNDARY 
SIDE1 , 2 
SIDE2 , 2 

★STEP,  NLGE0M,CYCLE=10 
★STATIC, PT0L=l.E-4, DIRECT 
1.  ,  1. 

★END  STEP 
★STEP , NLGEOM 

★STATIC, PT0L=l.E-4  , DIRECT 
1.  ,1. 

★MODEL  CHANGE,  REMOVE 
MANDREL, SFIT 
★END  STEP 


fa)  interference-fit  (b)  separation 

Figure  3.  A  two-dimensional  finite  element  model. 
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(  r-  cw  /  c  b  -  a  j 


Figure  4.  The  effect  of  interference  on  the  hoop 
stresses  using  rigid  mandrels. 


kj  the  bore. 


TECHNICAL  REPORT  INTERNAL  DISTRIBUTION  LIST 

NO.  OF 
COPIES 


CHIEF,  DEVELOPMENT  ENGINEERING  BRANCH 

ATTN:  SMCAR-CCB-D  1 

-DA  1 

-DC  1 

-DM  1 

-DP  1 

-DR  1 

-DS  (SYSTEMS)  1 

CHIEF,  ENGINEERING  SUPPORT  BRANCH 

ATTN:  SMCAR-CCB-S  1 

-SE  1 

CHIEF,  RESEARCH  BRANCH 

ATTN:  SMCAR-CCB-R  2 

-RA  1 

-RM  1 

-RP  1 

-RT  1 

TECHNICAL  LIBRARY  5 

ATTN:  SMCAR-CC8-TL 

TECHNICAL  PUBLICATIONS  &  EDITING  UNIT  3 

ATTN:  SMCAR-CCB-TL 

DIRECTOR,  OPERATIONS  DIRECTORATE  1 

ATTN:  SMCWV-00 

DIRECTOR,  PROCUREMENT  DIRECTORATE  1 

ATTN:  SMCWV-PP 

OIRECTOR,  PRODUCT  ASSURANCE  DIRECTORATE  1 

ATTN :  SMCWV-QA 


i 

!  NOTE:  PLEASE  NOTIFY  DIRECTOR,  BENET  LABORATORIES,  ATTN:  SMCAR-CCB-TL,  OF 

[  ANY  ADORESS  CHANGES. 


TECHNICAL  REPORT  EXTERNAL  DISTRIBUTION  LIST 


NO.  OF 
COPIES 

ASST  SEC  OF  THE  ARMY 

RESEARCH  AND  DEVELOPMENT 

ATTN:  DEPT  FOR  SCI  AND  TECH  1 

THE  PENTAGON 

WASHINGTON,  O.C.  20310-0103 
ADMINISTRATOR 

DEFENSE  TECHNICAL  INFO  CENTER 
ATTN:  DTIC-FDAC  12 

CAMERON  STATION 
ALEXANDRIA,  VA  22304-6145 


COMMANDER 
US  ARMY  ARDEC 

ATTN:  SMCAR-AEE  1 

SMCAR-AES ,  BLDG.  321  1 

SMCAR-AET-0 ,  BLDG.  351N  1 

SMCAR-CC  1 

SMCAR-CCP-A  1 

SMCAR-FSA  1 

SMCAR-FSM-E  1 

SMCAR-FSS-D ,  BLDG.  94  1 


SMCAR-IMI-I  (STINFO)  BLDG.  59  2 

PICATINNY  ARSENAL,  NJ  07806-5000 

DIRECTOR 

US  ARMY  BALLISTIC  RESEARCH  LABORATORY 
ATTN:  SLC8R-DD-T ,  BLDG.  305  1 

ABERDEEN  PROVING  GROUND,  MD  21005-5066 

DIRECTOR 

US  ARMY  MATERIEL  SYSTEMS  ANALYSIS  ACTV 
ATTN:  AMXSY-MP  1 

ABERDEEN  PROVING  GROUND,  MD  21005-5071 

COMMANDER 
HQ,  AMCCOM 

ATTN:  AMSMC-IMP-L  1 

ROCK  ISLAND,  IL  61299-6000 


NO.  OF 
COPIES 

COMMANDER 

ROCK  ISLAND  ARSENAL 

ATTN:  SMCRI-ENM  1 

ROCK  ISLAND,  IL  61299-5000 

DIRECTOR 

US  ARMY  INDUSTRIAL  BASE  ENGR  ACTV 
ATTN:  AMXIB-P  1 

ROCK  ISLAND,  IL  61299-7260 

COMMANDER 

US  ARMY  TANK-AUTMV  R&D  COMMAND 
ATTN:  AMSTA-DDL  (TECH  LIB)  1 

WARREN,  MI  48397-5000 

COMMANDER 

US  MILITARY  ACADEMY  1 

ATTN:  DEPARTMENT  OF  MECHANICS 
WEST  POINT,  NY  10996-1792 

US  ARMY  MISSILE  COMMAND 
REDSTONE  SCIENTIFIC  INFO  CTR  2 

ATTN:  DOCUMENTS  SECT,  BLDG.  4484 
REDSTONE  ARSENAL,  AL  35898-5241 

COMMANDER 

US  ARMY  FGN  SCIENCE  AND  TECH  CTR 
ATTN :  DRXST-SD  1 

220  7TH  STREET,  N.E. 

CHARLOTTESVILLE,  VA  22901 

COMMANDER 

US  ARMY  LABCOM 

MATERIALS  TECHNOLOGY  LAB 

ATTN:  SLCMT-IML  (TECH  LIB)  2 

WATERTOWN,  MA  02172-0001 


NOTE: 


PLEASE  NOTIFY  COMMANDER,  ARMAMENT  RESEARCH,  DEVELOPMENT,  AND  ENGINEERING 
CENTER,  US  ARMY  AMCCOM,  ATTN:  BENET  LABORATORIES,  SMCAR-CC8-TL , 
WATERVLIET,  NY  12189-4050,  OF  ANY  ADDRESS  CHANGES. 
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NO.  OF 
COPIES 

COMMANDER 

US  ARMY  LABCOM,  ISA 

ATTN:  SLCIS-IM-TL  1 

2800  POWDER  MILL  ROAD 
ADELPHI,  MD  20783-1145 

COMMANDER 

US  ARMY  RESEARCH  OFFICE 

ATTN:  CHIEF,  IPO  1 

P.O.  BOX  12211 

RESEARCH  TRIANGLE  PARK,  NC  27709-2211 
DIRECTOR 

US  NAVAL  RESEARCH  LAB 

ATTN:  MATERIALS  SCI  &  TECH  DIVISION  1 

CODE  26-27  (DOC  LIB)  1 

WASHINGTON,  D.C.  20375 


NO.  OF 
COPIES 

COMMANDER 

AIR  FORCE  ARMAMENT  LABORATORY 
ATTN:  AFATL/MN  1 

EGLIN  AFB,  FL  32542-5434 

COMMANDER 

AIR  FORCE  ARMAMENT  LABORATORY 
ATTN:  AFATL/MNF 

EGLIN  AFB,  FL  32542-5434  1 

METALS  AND  CERAMICS  INFO  CTR 
BATTELLE  COLUMBUS  DIVISION 
505  KING  AVENUE 

COLUMBUS,  OH  43201-2693  1 


NOTE: 


PLEASE  NOTIFY  COMMANDER,  ARMAMENT  RESEARCH,  DEVELOPMENT,  AND  ENGINEERING 
CENTER,  US  ARMY  AMCCOM,  ATTN:  BENET  LABORATORIES,  SMCAR-CCB-TL , 
WATERVLIET,  NY  12189-4050,  OF  ANY  ADDRESS  CHANGES. 


